AEDC-TR-94-14 


Analysis  of  Planar  Laser-Induced  Fluorescence 
Images  Obtained  During  Shakedown  Testing  of 
the  AEDC  Impulse  Facility 

W.  M.  Ruyten,  W.  D.  Williams,  and  F.  L.  Heltsley 
Micro  Craft  Technology/AEDC  Operations 

March  1995 

Final  Report  for  Period  October  1 993  —  September  1 994 


Approwd  tor  pidilie  ratasM;  datrlbutian  la  unNmltad. 


ARNOLD  ENGINEERING  DEVELOPMENT  CENTER 
ARNOLD  AIR  FORCE  BASE,  TENNESSEE 
AIR  FORCE  MATERIEL  COMMAND 
UNITED  STATES  AIR  FORCE 


NOTICES 


When  U.  S.  Government  drawinp,  specifications,  or  other  data  are  used  for  any  purpose  other  than  a 
definitely  related  Government  procurement  operation,  the  Government  thereby  incurs  no  responsibili^ 
nor  any  obligation  whatsoever,  and  the  fact  that  the  Government  may  have  formulated,  furnished,  or  in 
any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be  regarded  by  implication  or 
otherwise,  or  in  any  manner  licensing  the  holder  or  any  other  person  or  coiporation,  or  conveying  any 
rights  or  permission  to  manufacture,  use,  or  sell  any  patented  invention  that  may  in  any  way  be  related 
thereto. 


Qualified  users  may  obtain  copies  of  this  report  from  the  Eiefense  Technical  Information  Center. 


References  to  named  commercial  products  in  this  report  ate  not  to  be  considered  in  any  sense  as  an 
endorsement  of  the  product  by  the  United  States  Air  Force  or  the  Government. 


This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (PA)  and  is  releasable  to  the  National 
Technical  Information  Service  (NTIS).  At  NTIS,  it  will  be  available  to  the  general  public,  including 
foreign  nations. 


APPROVAL  STATEMENT 


This  report  has  been  reviewed  and  approved. 


Technical  Director 

Space  and  Missile  Systems  Tbst  Division 


Approved  for  publication; 


FOR  THE  COMMANDER 


CONRAD  M.  RICHEY,  LtCol,  USAF 
Chief,  Space  and  Missile  Test  Division 
Operations  Directorate 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMBNo.  0704-0188 

riporbng  burden  for  Ihit  oaleclnn  of  InTunnKion  le  ertimeted  lo  average  1  Nmt  par  retponaa.  eicludlrig  the  Ibne  for  reulaiMing  InatructiorTfr,  aNr^tbig  axiating  dela  eoureea. 
gathering  end  maintatring  the  data  needed,  and  compleitng  and  itviavnng  the  collacclDn  d1  Infaimatlon  Send  cornmenta  ragercfeno  ihia  biadan  esianaie  ar  any  other  aipect  of  thit 
eofeKtion  of  nfarmauan,  indudng  wggeeiiara  lor  reducing  ihia  bitfdan.  Co  WaBhlngten  Haadquiriere  SerwIceB.  Dtrectorata  for  Infomtalkon  Oparaeona  and  Heporta,  1218  JelferBOfi 
Devil  Hqhway.Siiie  12M,  Arlmgtcn.  VA  22202'43O2,  ind  lotheOfflca  of  Monaganwnl  and  Budgat*  Paporwork  Haductlon  Pro^t  (07044>188},  Washington,  DC  2060S. 

1  AGENCY  USE  ONLY  (Leova  blank)  2  REPORT  DATE  3  REPORT  TYPE  AND  OATES  COVERED 

March  1995  Final  -  October  1993  -  September  1994 

4.  TITIE  AND  SUBTITlf 

Analysis  of  Planar  Laser-Induced  Fluorescence  Images  Obtained  During  Shakedown 
Testing  of  the  AEDC  Impulse  Facility 

S.  FUNDING  NUMBERS 

0129 

8.  AUTHOR(E> 

W.  M.  Ruyten,  W.  D.  Williams,  and  F.  L.  Heltsley 

Micro  Craft  Technology,  AEDC  Operations 

7  rERFOnUIMU  □neANIZATION  HAMEISI  AMO  ADOfl ESSIES) 

Arnold  Engineering  Development  Center/ DOS 

Air  Force  Materiel  Command 

Arnold  Air  Force  Base,  TN  373S9-4000 

e.  PERFORMING  ORGANIZATION 

IREPORT  NUMBH) 

AEDC-TR-94-14 

9  SPONSOSMOMONnORING  AGENCY  NAMESISI  AND  AODRESSIES) 

Arnold  Engineering  Development  Center/DOS 

Air  Force  Materiel  Command 

Arnold  Air  Force  Base,  TN  37389-4000 

10  SPONSORWO/MONITORINO 

AGENCY  REPORT  NUMBER 

11  SUPPLEMENTARY  NOTES 

Available  in  Defense  Technical  Information  Center  (DTIC). 

12A  DISTmaUTfON/AVAILABILfTY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

12D  DISTRIBUTION  CODE 

13  ABSTRACT  (Maxanura  200  wnnlal 

Following  a  proof-of-princLple  demonstration  of  Planar  Laser-Induced  Ruorescence  (PUF)  in  the  AEDC  Impulse  Facility 
in  FY93,  a  set  of  four  PLIF  images  was  obtained  during  FY94  runs  of  the  facility.  The  images  were  obtained  aw^  from 
the  nose  region  of  a  spherically  blunt  cone,  using  excitation  of  nitric  oxide  in  the  flow.  Two  laser  sheet  geometries  were 
employed  and  two  excitation  wavelengths  were  used.  A  procedure  was  developed  whereby  calibration  images  allow  the 
remapping  of  the  PLIF  images  to  test  article  coordinates.  Using  this  technique,  it  was  shown  that  the  bow  shock  positions 
determined  by  PLIF  agree  well  with  the  shock  positions  measured  using  schlieren  photography.  More  importantly,  the  set 
of  PLIF  images  constitutes  an  initial  database  against  which  Computational  Fluid  Dynamics  calculations  of  high  enthalpy 
flows  can  be  compared.  Initial  results  of  such  a  comparison,  obtained  through  Conqrutational  Flow  Imaging,  are  reponed. 
Recommendations  are  given  for  future  PUF  ^plications  at  the  Impulse  Facility. 

14  SUBJECT  TU1MS 

flow  imaging  Computational  Flow  Imaging 

planar  laser-induced  fluorescence  imaging  processing 

hypersonic  flow  nitric  oxide 

IS.  NUMBER  OF  PAGES 

64 

10  PRICE  CODE 

17  SECURITY  CLASariCATlON  18  SECURfTY  CLA88FICATION  19  SECURITY  CLASSIFICATION 

OF  REPORT  OF  THB  PAGE  OF  ASSTRACT 

UNCLASSIFIED  UNCLASSIFIED  UNCLASSIFIED 

20  LMITATION  OF  ABSTRACT 

SAME  AS  REPORT 

,  Elindvd  Form  2M  (Rflv.  S-ftfi) 

COMPUTiR  GENERATED  PreKribad  by  ANSI  Scd  ZSSlB 


2M-102 


AEDC-TR-94-14 


PREFACE 

The  work  reported  here  was  done  at  the  Arnold  Engineering  Development  Center  (AEDC), 
Air  Force  Materiel  Command  (AFMC),  Arnold  AFB,  TN.  TTie  AF  project  manager  was  Gary  R. 
Mattasits;  the  manuscript  was  submitted  for  publication  on  March  1, 199S 

The  authors  acknowledge  the  contributions  of  R.  L.  Clippard,  J.  H.  Jones.  B.  J.  McClure,  J.  R. 
Maus,  L.  L.  Price,  S.  G.  Rock,  M.  S.  Smith,  and  J.  T.  Tarvin  in  various  stages  of  this  work. 
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1.0  INTRODUCTION 

Planar  Laser-Induced  Fluorescence  (PLIF)  imaging  has,  in  recent  years,  been  shown  to  be  a 
powerful  non-intrusive  diagnostic  technique  for  the  study  of  aerodynamic  flow  phenomena.  In 
particular,  PLIF  has  the  potential  of  rendering  instantaneous  measurements  of  temperature  and  den¬ 
sity  fields  in  the  flow.  However,  the  application  of  the  technique  in  a  complex  facility  poses  many 
challenges,  from  the  perspectives  of  experimental  complexity  as  well  as  data  analysis.  Thus  it  was 
a  signiHcant  achievement  when  it  was  demonstrated  in  FY93  that  it  is  possible  to  obtain  PI  TP  im¬ 
ages  in  the  AEDC  Impulse  Shock  Tunnel.  In  particular,  PLIF  images  were  obtained  in  the  nose 
region  of  a  blunt  cone  test  article. 

Here  the  development  of  image  analysis  techniques  for  the  quantitative  interpretation  of 
such  PLIF  images  is  reported.  In  particular,  the  analyses  of  the  PLIF  images  which  were  obtained 
during  FY94  shakedown  testing  and  evaluation  of  the  Impulse  Facility  (Runs  14, 15, 16,  and  19) 
are  described.  These  images  were  also  obtained  with  a  blunt  cone  test  article,  but  away  from  the 
nose  region. 

Two  distinct  categories  of  image  analysis  can  be  distinguished;  (1)  normalization  and  scaling 
and  (2)  derivation  of  flow  properties.  Analysis  in  the  first  category  does  not  require  specific 
knowledge  of  the  flow  under  study,  but  is  indispensable  for  the  quantitative  interpretation  of  the 
images.  It  involves  correcting  the  image  for  imperfections  in  the  camera  and  the  laser  sheet,  and  it 
involves  geometric  scaling  to  relate  image  pixel  coordinates  to  test  article  coordinates.  This  repent 
is  concerned  mostly  with  analysis  in  this  category. 

Work  in  the  second  category  has,  for  this  study,  taken  the  form  of  Computational  Flow 
Imaging  (CFI).  That  is,  synthetic  PUF  images  were  generated  from  Computational  Fluid 
Dynamics  (CFD)  solutions  of  the  Impulse  Facility  flow  around  the  test  article,  and  these  were 
compared  with  the  experintental  PLIF  images  after  normalization  and  scaling.  The  results  of  the 
CFI-PLIF  analysis  have  been  documented  in  Ref.  1  and  will  be  described  only  briefly  in  this  report. 

2.0  SUMMARY  OF  EXPERIMENTAL  RESULTS 

Figure  I  shows  a  schematic  of  the  experimental  setup  (Ref.  2).  Figure  2  shows  a  schematic 
of  the  test  article  used  in  the  PLIF  studies.  A  pulsed,  wavelength-tunable  UV  laser  beam  was 
formed  into  a  thin  sheet  and  directed  vertically  into  the  Impulse  Tunnel  through  an  optical  port  at 
the  top  of  the  facility.  The  width  of  the  sheet  at  the  location  of  the  test  article  was  approximately 
100  mm;  the  thickness  of  the  sheet  was  less  than  1  mm.  The  laser  was  tuned  to  a  particular  absorp¬ 
tion  line  of  the  nitric  oxide  (NO)  molecule.  Two  transitions  were  used  in  this  wofk:  (1)  the  225.716 
nm  Rii+<52i(13.5)  transition  of  the  Gamma  (0-0)  band  (Runs  14, 15,  an  16);  and  (2)  the  225.987 
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nm  R]  i+Q2|(7.S)  transition  of  the  same  band  (Runs  17, 18,  and  19).  These  lines  were  chosen  be¬ 
cause  they  provide  relatively  strong  fluorescence  signals  and  because  they  are  well  isolated  spec¬ 
trally  from  other  absorption  features  in  the  Gamma  (0,0)  band. 

The  resulting  fluorescence  signals  were  recorded  using  an  image-intensified  CCD  camera, 
with  an  array  size  of  578  x  384  pixels  and  a  14-bit  digitizer  readout.  The  camera  was  mounted  to 
the  side  of  the  tuimel,  with  a  line  of  sight  at  an  angle  of  about  24  deg  off-perpendicular  to  the 
centerline  of  the  flow.  Perpendicular  viewing  of  the  laser  sheet  would  have  been  preferable,  but  the 
optical  ports  perpendicular  to  the  flow  were  used  by  a  schlieren  imaging  system.  UV-pass  filters 
were  mounted  on  the  camera  to  minimize  the  unintentional  imaging  of  scattered  laser  light  and 
spontaneous  (non-laser-induced)  emissions  from  the  flow.  Moreover,  a  short  gate  width  (25  nsec 
for  Run  1 4;  SO  nsec  for  all  other  images)  was  chosen  for  the  intensifier  on  the  CCD  camera  to  limit 
the  contribution  of  spontaneous  flow  emission.  This  gate  time  is  also  short  enough  that  the  resulting 
image  is  essentially  frozen  in  time:  at  a  typical  flow  velocity  of  4500  m/sec,  a  50  nsec  gate  would 
cause  maximum  blurring  of  0.23  mm;  the  actual  blurring  should  be  even  less  as  a  result  of  centering 
of  the  laser  pulse  within  the  gate  and  collisional  quenching  of  the  laser-excited  state. 

Despite  the  short  gate  width  and  the  spectral  Altering  of  the  signal,  spontaneous  flow  emis¬ 
sion  was  found  to  dominate  the  signal  in  the  nose  region  of  the  test  article  in  an  initial  PLIF  survey 
for  Run  13  (see  Ref.  2,  Fig.  1).  Thus,  emphasis  in  this  work  was  placed  on  obtaining  high-quality 
PLIF  images  a  few  inches  downstream  from  the  nose  region. 

Six  Impulse  Facility  runs  were  conducted  during  this  study  (Runs  14  through  19).  The  PLIF 
system  was  operated  for  each  of  these  runs.  However,  no  triggers  were  received  for  runs  17  and  18 
so  that  no  PLEP  images  are  available  for  these.  In  Runs  14  and  IS,  the  laser  sheet  was  aligned 
parallel  to  the  model  surface,  so  as  not  to  be  intercepted  by  the  test  article  (Fig.  3a).  For  Runs  16 
and  19,  the  laser  sheet  was  aligned  along  the  meridian  plane  of  the  test  article,  impinging  on  the 
conical  surface  of  the  test  article  as  shown  in  Fig.  3b. 

Analysis  of  the  PLIF  images  for  Runs  14  and  15  is  complicated  by  the  fact  that  it  appears 
that,  in  both  images,  the  position  of  the  laser  sheet  moved  significantly  from  the  original  alignment 
position.  After  noting  this  problem  on  Run  IS,  the  laser  sheet  optics  were  fastened  more  securely 
for  subsequent  runs,  and  no  changes  in  alignment  were  seen  as  a  result  of  the  Impulse  Facility  fir¬ 
ings  for  Runs  16  and  19. 
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3.0  COMPARISON  WITH  SCHLIEREN  IMAGES 

Schlieren  images  were  obtuned  for  most  of  the  Impulse  Facility  runs.  Although  not 
specifically  intended  for  this  purpose,  they  were  used  to  verify  that  shock  locations  in  the  PLIF 
images  agree  with  those  in  the  schlieren  images  for  the  same  runs.  This  was  accomplished  as 
follows.  First,  a  schlieren  image  (the  one  for  Run  16  is  shown  in  Fig.  4)  was  put  on  a  digitizing 
tablet,  and  two  series  of  points  were  selected  manually:  one  along  the  model  surface  of  the  test 
article,  the  other  along  the  bow  shock. 

Next,  the  digitizer  tablet’s  model  surface  coordinates  were  least-squares  curve-fitted  to  the 
known  mathematical  definition  of  the  model  surface,  which  is  given  by: 


_)  ±*Ir^ -x^;  x<-f?sina, 

^  (  ±(i?  +  jcsina  )/cosa;  x^~I?sina, 

where  a  =  10°  is  the  half  angle  of  the  blunt  cone,  R  =  0.7  in.  is  the  radius  of  the  spherical  nose 
section,  and  the  origin  (x  =  y  =  0)  is  located  at  the  center  of  the  spherical  nose  section.  The  scaling 
transformation  from  test  article  coordinates  to  digitizer  tablet  coordinates  is  determined  by  four  fit 
parameters;  a  linear  zoom  factor,  a  horizontal  and  a  vertical  offset,  and  a  rotation  angle.  After 
obtaining  these  fit  parameters,  the  digitizer  tablet's  coordinates  for  the  shock  position  could  be 
transformed  to  test  article  coordinates  by  applying  the  inverse  transformation. 

Finally,  the  schlieren-derived  shock  positions  were  cast  into  parametric  form  by  fitting  the 
test  article  coordinates  of  the  shock  to  the  empirical  equation 

y  =  ±1^  i(a  +X  +  CX  )J  ,  (2) 

where  a.  b.  c.  and  X  are  the  parameters  of  the  least-squares  fit.  As  an  example.  Fig.  5  shows  the 
theoretical  model  surface,  the  fitted  shock  location,  and  the  rescaled  digitizer  coordinates  for  the 
schlieren  image  obtained  during  Run  16.  Fitted  values  of  the  parameters  a,  b,  c,  and  X  for  Runs  14, 
IS,  and  16  are  given  in  Table  I  (valid  for  x  and  y  in  mm).  These  parameters  are  highly  correlated. 
Thus  significant  differences  between  the  fit  parameters  are  not  necessarily  indicative  of  large  vari¬ 
ations  in  shock  position.  This  is  evident  in  Fig.  6,  in  which  the  fitted  schlieren  shock  positions  are 
shown  for  Runs  14,  IS,  and  16. 
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Finally  it  is  noted  that  the  above  equation  for  the  schlieren  shock  position  is  symmetric 
around  the  centerline  (y  =  0)  by  definition.  However,  slight  deviations  from  symmetry  were  seen 
in  the  digitized  shock  positions  (Fig.  3)  as  well  as  in  the  schlieren  images  themselves  (  Fig.  4).  Thus 
it  appears  that  some  flow  asymmetry  was  present. 

4.0  TRANSFORMATION  FROM  CAMERA  TO  TEST  ARTICLE  COORDINATES 

To  establish  the  coordinate  transformation  between  camera  pixel  coordinates  and  test  article 
coordinates,  two  calibration  images  were  typically  obtained  prior  to  an  Impulse  Facility  run.  The 
Srst  was  an  image  of  a  calibration  target,  consisting  primarily  of  a  rectangular  array  of  marks;  the 
second  image  was  either  a  front-lit  or  back-lit  image  of  the  actual  test  article  in  which  the  model 
edges  can  be  clearly  discerned. 

The  calibration  targets  were  placed  at  the  anticipated  positions  of  the  laser  sheets.  In  this 
work,  a  target  consisting  of  an  array  of  equally  spaced  rectangular  dots  was  used.  Some  horizontal 
and  vertical  crosshair  lines  were  marked  on  the  target  as  well.  Figures  7  and  8  show  examples  of 
the  resulting  target  images  for  Runs  15  and  16.  In  Fig.  7,  the  target  was  placed  at  the  pre-run  laser 
sheet  position  shown  in  Fig.  2a;  that  is,  it  was  aligned  parallel  to  the  model  surface.  In  Fig.  8,  the 
target  was  cut  at  a  10-deg  angle  and  placed  on  top  of  the  test  article,  aligned  along  the  meridian 
plane  of  the  test  article  (Fig.  3b).  In  both  cases,  the  horizontal  alignment  of  the  target  relative  to  the 
test  article  was  achieved  by  centering  the  copper-stainless  oint  on  the  test  article  (3  in.  from  the 
model  nose)  between  two  of  the  vertical  cross  hairs  in  the  target  (in  the  case  of  Fig.  7,  this  alignment 
was  facilitated  by  a  small  hole  in  the  target).  Figure  9  shows  the  back-lit  image  of  the  test  article 
for  Run  15.  Because  the  model  edge  is  visible  in  Fig.  8,  no  separate  back-lit  image  was  required. 

Because  the  camera  is  mounted  about  24  deg  off-perpendicular,  the  horizontal  features  in  the 
calibration  images  of  Figs.  7  and  8  are  seen  to  converge  to  the  left.  This  results  in  a  nonlinear  trans¬ 
formation  between  ob  ect  and  image  coordinates.  To  describe  this  transformation,  the  analysis  was 
generalized  to  the  case  where  the  off-perpendicular  viewing  angle  of  the  camera  is  given  by  an 
arbitrary  angle,  Y,  and  the  viewing  direction  is  pitched  away  from  the  horizontal  plane  by  an  (small, 
unintended)  angle,  5  (Fig.  10).  Furthermore,  the  camera  was  assumed  to  be  pointed  (and  focused) 
at  some  reference  point  (JC(>,yo.Zo)’  ^  removed  by  a  distance,  D,  from  this  reference  point,  and  to 

have  associated  with  it  an  effective  zoom  factor,  p.  Neglecting  aberration,  the  image  of  an  arbitrary 
point  ix,y,z)  may  now  be  expressed  in  terms  of  camera  pixel  coordinates  (Xp.yp)  as 
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P  (w  ■  g,) 
D-  (it  ■  ej 

P  (w  •  gp 


(3) 


where  iy,  and  e^.  are  mutually  orthonormal  unit  vectors  given  by 

=  ( cosy, 0.  -  sin y) 
gy  =  (-siny  sin6,  cos5,  -cosy  sin6) 

e„  =  (siny  cos 6, sin 6, cosy  cosS)  (4) 

(these  vectors  are  oriented,  respectively,  along  the  x  and  y  axes  of  the  camera  and  along  the  camera 
viewing  direction),  and  w  is  the  vector  given  by 

it  =  (jfo-Jc.yo-y,Zo-z).  (5) 

The  parameters  which  must  be  determined  on  the  basis  of  the  calibration  images  are  the  angles  y 
and  5,  the  distance  A  the  zoom  factor  p,  and  the  components  of  the  reference  point,  xq,  y^,  and  zq- 
This  was  accomplished  through  a  least-squares  fitting  program,  which  takes  as  input  the  image 
centers  of  the  marks  in  the  calibration  target  and  the  image  coordinates  of  the  model  edges.  Both 
sets  of  coordinates  were  determined  from  the  calibration  images  through  a  semi-automatic  proce¬ 
dure  in  which  the  analyst  pulls  the  image  up  on  a  computer  screen  and  manually  selects  points  of 
interest  (either  dot  centers  or  edge  points),  which  are  then  optimized  automatically  by  dot-centering 
or  edge-centering  algorithms  in  the  analysis  software  (further  details  can  be  provided  on  request). 

As  a  check  on  this  procedure,  it  is  good  practice  to  compare  the  input  coordinates  and  fitted 
coordinates.  Typically,  an  average  fit  error  on  the  order  of  a  one  pixel  separation  was  found  to 
result.  As  an  example.  Fig.  1 1  shows  the  input  and  litted  coordinates  of  the  centers  of  the  reference 
marks  from  the  calibration  image  in  Fig.  7.  (In  Fig.  7,  there  is  also  a  slight  vertical  misalignment 
between  the  camera  and  the  calibration  target;  this  could  be  compensated  for  by  including  a  “roll” 
angle  tp  in  the  above  equations.  However,  this  inclusion  was  determined  not  to  have  a  significant 
effect  on  the  eventual  remapping  of  the  PLIF  images.)  Similarly,  Fig.  1 2  shows  the  projected  model 
position  (edges,  centerline,  and  model  joint  at  3  in.  from  the  nose)  as  an  overlay  on  Fig.  8.  Note  the 
good  alignment  of  the  projected  features  with  those  in  the  image. 
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5.0  REMAPPING  OF  A  PLIF  IMAGE  TO  TEST  ARTICLE  COORDINATES 


Having  established  the  transfonnation  from  test  article  coordinates  to  image  pixel 
coordinates,  two  approaches  can  be  taken.  The  simplest  of  these  is  to  pro  ect  a  known  feature  in 
the  flow  onto  the  PLIF  image.  This  approach  proved  very  useful  in  checking  for  misalignment 
effects  of  the  actual  PLIF  image  relative  to  the  pre-run  calibration  images.  It  was  also  instrumental 
in  establishing  the  actual  position  of  the  laser  sheet  for  Runs  14  and  15,  where  motion  of  the  laser 
sheet  during  the  Impulse  Facility  run  appears  to  have  taken  place  (see  Sec.  8.0). 


More  generally,  the  coordinate  transformation  between  test  article  coordinates  and  image 
coordinates  can  be  used  to  remap  a  PLIF  image  to  test  article  coordinates,  and  to  undo  the  nonlinear 
mapping  which  is  introduced  by  the  ofr-petpendicular  viewing  direction  of  the  camera. 
Specifically,  when  the  position  of  the  PLIF  laser  sheet  is  known  (for  example,  z  =  0  for  Runs  16  and 
19,  and  z  =  for  Runs  14  and  15,  where  a,  and  b,  are  known),  it  is  possible  to  invert  Eqs.  (3, 
4,  and  5),  yielding  test  article  coordinates  (jc,y,z)  in  terms  of  image  pixel  coordinates  (.Xp,yp). 
Schematically,  this  may  be  expressed  as 

ix.y.z)  =  T'  (i^,yp) .  (6) 


However,  to  remap  an  image,  the  inversion  expressed  in  Eq.  (6)  does  not  need  to  be  performed 
explicitly.  Instead,  a  new  image  grid  is  defined,  for  example,  along  the  meridian  plane  of  the  test 
article.  In  order  to  remap  the  PLIF  image  onto  this  newly  defined  image  grid,  the  virtual  “image” 
of  each  grid  cell  in  the  new  grid  is  calculated  using  Eqs.  (3,  4,  and  5).  The  result  is  a  set  of  four 
(fractional)  pixel  coordinates  in  the  original  PLIF  inuige  (Fig.  13),  which  may  span  less  than  one 
pixel  (Fig.  13a),  about  one  pixel  (Fig.  13b),  or  a  number  of  pixels  (Fig.  13c)  in  the  PLIF  image  (or 
one  or  more  of  the  "imaged"  comers  may  lie  outside  of  the  boundary  of  the  PLIF  image).  The 
remapped  PLIF  intensity  assigned  to  the  pixel  in  the  new  grid  is  now  given  by  a  weighted  average 
of  pixel  values  in  the  original  PLIF  image.  That  is. 


'rtmapptd 


=  1/, 


(f.;) 

original 


XF 


(U) 


(7) 


where  the  summation  is  over  all  pixels  (ij)  which  are  inside  the  imaged  regions  in  Fig.  13,  and 
are  the  enclosed  fractions  of  the  pixels  in  the  original  image.  It  may  be  noted  that  the  resulting 
remapped  image  may  have  about  the  same  spatial  resolution  as  the  original  image  (as  in  Hgs.  13a 
and  13b)  or  less  spatial  resoludon  (as  in  Fig.  13c).  Remapping  of  this  sort  has  been  applied  to  each 
of  the  PLIF  images  and  is  discussed  in  Secs.  7.0  and  8.0. 


10 


AEDC-TR-94-14 


6.0  BACKGROUND  SUBTRACTION  AND  FLAT-FIELD  CALIBRATION 

Two  corrections  which  were  performed  on  the  raw  PLIF  images  before  further  processing 
were  a  dark  count  subtraction  and  a  flat-field  correction.  The  dark  count  subtraction  is  necessary 
because  even  in  the  absence  of  an  optical  signal,  a  nonzero  count  rate  is  registered  (camera  dark 
count).  The  dark  count  rate  varies  by  only  a  few  counts  across  the  image  and  is  assumed  to  be 
given  by  a  constant  value.  This  value  can  easily  be  determined  from  an  image  with  the  camera 
lens  blocked. 

Next,  flat-rield  corrections  were  performed.  This  was  done  by  dividing  an  image  (after  dark 
count  subtraction)  by  a  flat-rield  calibration  (FFC)  image,  which  is  essentially  the  camera  response 
for  a  uniform  illumination  of  the  ICCD  array.  Figure  14  shows  the  FFC  image  which  was  used  for 
this  purpose.  This  calibration  was  performed  using  a  continuous  light  source  and  with  a  rather  long 
gate  time  on  the  intensifier  array  (Ref.  3).  Two  features  stand  out.  The  first  is  a  honeycomb  pattern 
spread  over  the  entire  image,  characterized  by  locally  reduced  sensitivities.  This  pattern  is  the  result 
of  the  hexagonal  structure  of  the  intensifier  array.  The  other  feature  is  a  more  global  variation  in 
camera  sensitivity  across  the  image.  For  the  FFC  file  displayed  in  Fig.  14,  this  global  variation 
amounts  to  a  variation  in  sensitivity  across  the  image  of  about  -25  to  25  percent. 

Subsequent  calibrations  of  the  camera  system  have  shown  that  the  FFC  response  depends,  to 
some  degree,  on  the  mode  of  operation  of  the  camera.  For  example,  Fig.  15  shows  a  typical  FFC 
image  obtained  with  a  short  gate  time  of  200  nsec  and  a  pulsed  laser  source.  This  FFC  response  is 
probably  more  typical  for  the  short  gate  width  operating  mode  than  the  FFC  response  in  Fig.  14, 
which  was  obtained  with  a  long  exposure  time  and  a  continuous  light  source.  Analysis  of  the  FFC 
calibration  data  is  still  in  progress.  In  all  cases  reported  here  where  an  FFC  was  performed,  the  old 
calibration  from  Fig.  14  was  used. 

Finally,  in  the  PLIF  images,  there  is  typically  a  small  background  signal,  spread  roughly 
uniformly  across  the  image.  This  background  was  presumed  to  be  the  result  of  diffuse  laser  scatter 
and/or  secondary  fluorescence  emission  (reabsorption  and  reemission  of  fluorescence  radiation  out¬ 
side  the  PLIF  laser  sheet).  Typically,  this  background  level  was  determined  by  calculating  the 
column-averaged  signal  outside  of  the  laser  sheet  and  performing  a  linear  tit  across  the  image.  For 
example,  for  the  PLIF  image  for  Run  14,  a  background  baseline  was  found  that  slopes  linearly  from 
6  counts  on  the  left  side  of  the  image  to  2  counts  on  the  right.  Following  dark  count  subtraction  and 
flat-field  correction,  such  baseline  signals  were  also  subtracted  from  the  images. 
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7.0  PLIF  IMAGES  FX)R  RUNS  16  AND  19 


7.1  RUN  16 

Because  they  ate  more  easily  inteipreted,  the  PLIF  images  for  Runs  16  and  19  are  discussed 
before  those  of  Runs  14  and  IS.  Figure  16  shows  the  raw  PLIF  image  obtained  for  Run  16.  The 
laser  sheet  was  directed  along  the  meridian  plane  of  the  test  article  (see  Bg.  2b).  Thus  there  is  no 
PLIF  signal  below  the  surface  of  the  model.  The  edges  of  the  laser  sheet  are  clearly  discernible. 
Figure  17  shows  the  same  image  in  which  the  model  surface,  the  model  oint  at  3  in.  from  the  nose, 
and  the  schlieren  shock  position  have  been  pro  ected.  This  shock  position  was  derermined  as 
described  in  Sec.  3.0.  Specincally,  the  dotted  line  in  Fig.  17  was  obtained  by  pro  ecting  the  curve- 
fitted  position  of  the  schlieren  shock  onto  the  PLIF  image.  Both  the  pro  ected  model  surface  and 
the  pro  ected  shock  match  the  corresponding  features  in  the  PLIF  image  well.  This  is  convincing 
proof  that  the  approach  to  obtaining  the  geometric  scaling  between  test  article  and  image 
coordinates  through  the  use  of  calibration  images  as  described  in  Sec.  4.0  is  sound.  In  particular, 
the  pre-run  scaling  transformation  could  be  applied  without  the  need  for  any  fine-tuning,  indicating 
that  none  of  the  optics  deviated  from  the  pre-run  alignment  during  the  Impulse  Facility  run. 

Clearly,  the  intensity  profile  of  the  laser  sheet  in  Figs.  15  and  16  is  nonuniform.  To  correct 
for  this,  the  average  PLIF  signal  was  calculated  in  the  free-stream  region  in  Fig.  16.  Essentially, 
this  profile  was  obtained  by  column  averaging  the  image  in  the  free-stream  region.  However,  the 
small  divergence  of  the  laser  sheet  was  accounted  for  (namely,  by  manually  selecting  the 
nonveitical  edges  of  the  laser  sheet  before  calculating  the  average  profile).  Also,  a  spline-smooth¬ 
ing  algorithm  was  employed  to  reduce  the  noise  on  the  resulting  profile,  which  is  shown  in  Fig.  18. 
The  horizontal  scale  in  Fig.  18  is  essentially  a  pixel  scale,  but  corrected  for  sheet  divergence.  The 
signal  outside  of  the  sheet  was  used  for  the  determination  of  the  background  signal. 

Following  this  procedure,  the  FLIP  image  was  energy-normalized  by  dividing  the  PLIF 
image,  roughly  column-wise  (while  accounting  for  the  slight  beam  divergence),  by  the  net  profile 
in  Fig.  18  (that  is,  after  subtracting  the  sloping  base  line).  Obviously,  this  normalization  is  only 
meaningful  where  the  laser  intensity  is  above  the  noise  level  of  the  image,  so  that  the  normalized 
image  has  to  be  truncated  at  the  laser  sheet  edges  (and,  of  course,  the  image  must  be  truncated  below 
the  model  surface,  since  there  is  no  PLIF  signal  below  the  edge  of  the  model). 

Finally,  after  performing  the  dark  count  subtraction,  the  flat-field  calibration,  the  back¬ 
ground  subtraction,  and  the  energy  normalization,  the  PLIF  image  was  remapped  to  test  article 
coordinates  using  the  procedure  described  in  Sec.  S.O.  The  result  is  shown  in  Fig.  19.  Also  shown 
in  Fig.  19  is  the  edge  of  the  test  article  and  the  model  oint  at  3  in.  from  the  nose. 
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Figure  20  shows  average  radial  profiles  of  the  PLIF  signal  through  the  center  of  Fig.  19  (at 
an  axial  position  of  83.0  mm),  as  well  as  average  profiles  to  the  left  and  to  the  right  of  the  center. 
Hie  profiles  were  obtained  by  averaging  the  PLIF  signal  over  28.5-mm-wide  bands,  each  spanning 
one-third  of  the  width  of  the  laser  sheet.  A  dimensionless  radial  coordinate  was  chosen  which  is 
zero  at  the  model  surface  and  one  at  the  shock  (as  determined  from  the  schlieren  image).  In  this 
way,  averaging  over  a  band  of  finite  width  does  not  cause  an  appreciable  blurring  of  the  shock 
position. 

Extensive  discussion  of  the  flow  features  in  Figs.  19  and  20  is  given  in  AIAA  paper  94-2621 
(Ref.  1).  In  particular,  the  comparison  between  this  processed  PLIF  image  and  the  synthetic  PLIF 
image  obtained  through  Computational  Flow  Imaging  is  discussed  in  Ref.  1.  Briefly,  there  is 
reasonable  agreement  between  the  experimental  and  theoretical  radial  profiles  around  the  shock 
position;  both  are  characterized  by  an  increase  in  PLIF  signal  when  passing  from  the  fiee-stream 
into  the  shock,  followed  by  a  subsequent  drop  in  signal  when  moving  even  closer  to  the  model 
surface.  However,  whereas  the  CFD/CFI  calculation  predicts  that  the  signal  should  continue  to 
decrease  toward  the  model  surface,  the  experimental  signal  rises  again  near  the  model  surface. 

This  discrepancy  has  not  been  resolved  and  may  be  explained  in  several  ways;  (1)  the  NO 
density  predicted  by  the  CFD  code  is  too  low  near  the  surface;  (2)  the  temperature  predicted  by  the 
CFD  code  is  too  high  near  the  surface;  or  (3)  the  experimental  PLIF  signal  close  to  the  model 
surface  contains  a  laser-induced  fluorescence  component  from  an  unidentified  flow  species  or  an 
unidentified  spectral  line  of  the  NO  molecule.  The  latter  possibility  seems  remote,  because  no  other 
tabulated  spectral  features  could  be  identified  at  the  excitation  wavelength  used  for  Runs  14,  15, 
and  16  (225.716  nm). 

(Following  the  presentation  of  Ref.  1,  further  CFD/CFI  simulations  have  been  performed. 
These  indicate  that,  unlike  assumed  ad-hoc  in  Ref.  1,  thermal  nonequilibrium  in  the  free-stream 
region  is  not  responsible  for  the  decrease  in  PUF  signal  when  passing  from  the  free  stream  into  the 
shock.  Thermal  nonequilibrium  (unequal  rotational  and  vibrational  temperatures  of  the  NO 
absorbers)  does  appear  to  be  present,  but  does  not  seem  to  significantly  affect  the  PLIF  signal  in  the 
transition  between  the  free  stream  and  the  shock.  However,  attenuation  of  the  laser  sheet  as  it 
propagates  through  the  free  stream  was  not  included  in  the  analysis  discussed  in  Ref.  1  and  has  now 
been  shown  to  cause  this  very  decrease  in  PLIF  signal  when  passing  from  the  free  stream  into  the 
shock.  But,  even  after  incorporating  the  effect  of  laser  attenuation  in  the  CFD/CFI  analysis,  the 
experimentally  observed  rise  in  PLIF  signal  when  approaching  the  model  surface  is  not  explained 
by  the  CFD/CFI  calculations;  a  more  complete  analysis  will  be  presented  at  a  liUer  time.) 


13 


AEDC-TR-94-14 


7.2  RUN  19 

To  reduce  speculation  that  an  accidental  spectral  coincidence  with  another  species  was  the 
reason  for  the  unexplained  signal  close  to  the  model  surface  in  Run  16,  the  laser  was  operated  at  a 
different  wavelength  (225.987  nra)  for  Runs  17,  18,  and  19.  Unfortunately,  no  triggers  were 
received  for  Runs  17  and  18.  The  PLIF  image  resulting  for  Run  19  is  shown  in  Engs.  21  and  22. 
Figure  21  shows  the  raw  image;  Fig.  22  shows  the  fully  processed  image,  which  may  be  compared 
directly  with  that  for  Run  16  in  Fig.  19. 

Population  fractions  of  the  laser-excited  state  in  NO,  absorption  cross  sections,  and 
fluorescence  yields  should  be  roughly  similar  for  the  225.987-nm  transition  and  that  used  for  the 
other  PLIF  images  (225.716  ora).  This  is  shown  in  Fig.  23,  which  shows  the  temperature 
dependence  of  the  PLIF  signals  at  the  two  wavelengths  assuming  equal  NO  number  densities  (the 
quantity  plotted  in  Fig.  23  is  the  product  of  the  NO  population  fraction  and  the  Einstein  B 
coefficient  for  the  respective  transitions,  assuming  equal  rotational  and  vibrational  temperatures). 
Also,  the  average  laser  pulse  energies  were  comparable  for  the  two  images.  In  fact,  the  camera  had 
been  moved  closer  to  the  tunnel  window,  which  should  have  resulted  in  a  gain  in  signal  of  about  a 
factor  2.5  for  Ron  19  compared  with  Run  16.  Nevertheless,  the  signal  level  (and  hence  the  signal- 
to-noise  ratio)  for  Run  19  is  much  lower  than  that  for  Run  16,  for  unknown  reasons. 

Still  it  appears  that  the  character  of  the  image  for  Run  1 9  is  very  similar  to  that  for  Run  16.  In 
paiticulai,  a  dark  band  can  be  seen  to  run  through  the  images  in  Fig.  22,  much  like  the  one  which 
runs  through  the  processed  PLIF  image  in  Fig.  19.  There  is  also  a  somewhat  smaller  dark  band  ust 
above  the  model  surface,  immediately  to  the  right  of  the  model  oint. 

Figure  24  shows  average  radial  profiles  at  the  center  of  the  PLIF  image  from  Eng.  22  (at  an 
axial  location  of  76.5  mm),  as  well  as  average  profiles  to  the  left  and  to  the  right  of  the  center.  These 
profiles  were  obtained  in  similar  fashion  to  those  shown  for  Run  16  in  Fig.  20,  namely  by  averaging 
over  three  vertical  regions,  each  with  a  width  equal  to  one-third  the  width  of  the  laser  sheet  (28.1 
mm).  The  shapes  of  the  profiles  in  Fig.  24  are  quite  similar  to  those  in  Fig.  22.  The  main  difference 
is  that  the  ratio  between  the  peak  signal  inside  the  shock  and  the  free-stream  value  signal  is  smaller 
(about  1 .4)  than  it  is  for  Run  16  (about  2.2).  This  trend  is  consistent  with  the  different  temperature- 
dependences  of  the  two  transitions  shown  in  Fig.  23.  For  a  more  thorough  quantitative  check,  a 
CFD  solution  for  Run  19  would  be  required. 
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8.0  PLIF  IMAGES  FOR  RUNS  14  AND  15 


8.1  RUM4 

FLIP  images  were  obtained  for  Runs  14  and  IS.  In  both  cases,  the  laser  sheet  was  aligned 
so  as  to  pass  in  front  of  the  model  surface  without  being  intercepted  by  the  test  article  (see  Fig.  3a). 
This  alignment  was  performed  so  that  the  laser  sheet  was  parallel  to  the  model  surface  and  separated 
horn  it  by  about  5  mm.  However,  it  appears  that,  for  both  Run  14  and  Run  15.  the  laser  sheet  moved 
during  the  Impulse  Facility  run,  as  argued  below. 

Figure  25  shows  the  raw  PLIF  image  for  Run  14.  The  overall  signal  level  is  rather  low. 
Some  natural  flow  emission  is  observed  on  the  left  side  of  the  image,  ust  downstream  from  the  nose 
cone.  Also,  with  proper  pseudo-color  rendering,  an  even  smaller  level  of  background  emission  can 
be  discerned  on  the  far  right-hand  side,  along  the  edges  of  the  test  article. 

Figure  26  shows  the  average  horizontal  intensity  profile  obtained  from  Fig.  25.  This  profile 
was  obtained  by  column  averaging  over  the  entire  height  of  the  signal.  Based  on  Fig,  26  a  linearly 
sloping  baseline  was  assigned  for  the  background  signal  (shown  as  the  dashed  line  in  Fig.  26),  and 
the  raw  image  was  normalized  by  subtracting  this  linearly  sloping  baseline  and  dividing  the  result¬ 
ing  image  by  the  average  profile  of  Fig.  25.  The  resulting  energy-normalized  image  is  shown  in 
Fig.  27. 

Several  observations  are  noteworthy.  First,  although  the  shape  of  the  high-intensity  signal 
region  appears  to  be  roughly  symmetric  around  the  centerline,  the  magnitude  of  the  signal  is  not. 
In  particular,  the  average  signal  in  the  boOom  half  of  the  image  is  only  about  60  percent  of  that  in 
the  top  half.  This  can  also  be  seen  in  Fig.  28,  which  shows  the  average  vertical  profile,  obtained 
simply  by  row  averaging  across  the  whole  width  of  the  energy-normalized  image  shown  in  Fig.  27, 
The  observed  asymmetry  is  consistent  with  the  hypothesis  that  the  laser  sheet  is  being  attenuated  as 
a  result  of  absorption  by  NO.  Work  is  in  progress  to  estimate  the  NO  number  density  on  the  basis 
of  this  observed  attenuation. 

The  second  observation  is  that  the  vertical  extent  of  the  high-intensity  signal  region  in  Figs. 
25  and  27  is  not  consistent  with  the  expected  size  of  this  region  based  on  the  original  alignment  of 
the  laser  sheet.  To  illustrate  this.  Fig.  29  shows  the  pro  ected  edges  of  the  test  article  (as  they  would 
have  been  seen  by  the  camera),  the  outline  of  the  shock  (as  determined  from  the  schlieren  image  for 
Shot  14),  and,  most  importantly,  the  pro  ection  of  the  intersection  of  the  laser  sheet  with  the  shock. 
This  pro  ection  was  calculated  by  assuming  that  the  shock  is  axially  symmetric  and  that  the  laser 
sheet  was  in  the  original  alignment  position  (parallel  to  the  model  surface  at  a  distance  of  S  mm). 
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On  the  basis  of  the  PLIF  images  from  Runs  16  and  19,  it  would  be  expected  that  the 
intersection  of  the  laser  sheet  with  the  shock  would  show  up  clearly  in  the  PLIF  image  (Rgs.  19  and 
22).  However,  the  height  of  the  predicted  intersection  in  Fig.  29  is  clearly  much  taller  than  the 
height  of  the  high-signal  region  in  Figs.  25  and  27.  The  most  plausible  assumption  for  this 
discr^ancy  appears  to  be  that  the  laser  sheet  moved  hrom  its  original  alignment  position  during  the 
Impulse  Facility  firing. 

To  test  this  assumption,  the  intersections  of  the  laser  sheet  and  the  shock  were  calculated  for 
a  number  of  sheet  positions.  For  example.  Fig.  29  shows  the  intersections  for  a  series  of  separations 
“d"  (measured  radially)  between  the  laser  sheet  and  the  model  surface,  while  keeping  the  sheet 
parallel  to  the  surface  (the  intersection  marked  by  a  solid  line  in  Fig.  29  corresponds  approximately 
to  the  one  which  would  have  resulted  for  the  original  alignment  of  the  laser  sheet).  To  establish  the 
optimum  fit  of  the  observed  intersection  to  the  measured  intersection,  a  least-squares  fit  was 
performed  on  hand-picked  paints  from  the  PLIF  image,  selected  along  the  top  and  bottom  edges  of 
the  high-signal  region  in  the  PLIF  image.  These  points,  along  with  the  fitted  intersection,  are  shown 
as  overlays  on  the  normalized  PLIF  image  in  Hg.  30.  The  laser  sheet  position  which  corresponds 
to  this  least-squares  fit  result  is  shown  in  Fig.  31,  along  with  the  model  surface,  the  shock,  and  the 
original  position  of  the  laser  sheet.  Note  that  this  sheet  position  has  moved  considerably  farther 
away  from  the  model  surface,  and  is  no  longer  parallel  to  the  surface. 

Having  established,  in  this  manner,  the  actual  position  of  the  laser  sheet,  it  was  possible  to 
remap  the  PLIF  image  to  a  meridian  plane  through  the  test  article.  The  result  is  shown  in  Fig.  32, 
both  for  the  top  half  of  the  PLIF  image  (Fig.  32a)  and  for  the  bottom  half  (Fig.  32b;  the  signal  in 
this  image  was  divided  by  0.6  to  compensate  approximately  for  the  lower  signal  as  a  result  of 
attenuation  of  the  laser  sheet).  Also  shown  in  Fig.  32  are  the  pro  ected  positions  of  the  shock,  the 
outline  of  the  model  surface,  and  the  vertical  oint  in  the  test  article  (Figs.  19  and  22). 

A  notable  difference  between  Figs.  32  and  22  is  that  the  ratio  between  the  PLIF  signal  inside 
the  shock  to  that  in  the  free  stream  is  much  larger  for  Run  14  than  for  Run  16,  where  it  is 
approximately  2.2.  Because  of  the  extremely  low-signal  level  in  the  low-intensity  region 
(presumed  to  be  free  stream)  for  Figs.  27  and  32,  an  accurate  value  of  the  “free-stream” 
fluorescence  cannot  be  obtained.  However,  the  ratio  is  certainly  much  larger  than  2.2.  Based  on 
Fig.  28,  this  ratio  is  stimated  to  be  about  20  (±5).  No  good  explanation  is  presently  available  for 
this  difference. 

RUN  15 

The  results  for  Run  15  are  very  similar,  except  that  a  much  stronger  signal  was  measured. 
Figure  33  shows  the  raw  PLIF  image.  Again,  there  is  a  high-intensity  region  with  a  shape 
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inconsistent  with  the  expected  intersection  of  the  laser  sheet  and  the  shock.  More  structure  is 
evident  in  Fig.  33  than  in  Figs.  25  and  27.  Also,  as  in  Fig.  25,  the  signal  level  in  the  bottom  part  of 
the  figure  is  lower  than  that  in  the  top  part,  leading  again  to  the  hypothesis  that  laser  attenuation  was 
present  as  a  result  of  absorption  by  NO. 

Figure  34  shows  the  average  horizontal  profile  (obtained  by  averaging  across  the  entire 
height  of  the  image),  and  Fig.  35  shows  the  energy-normalized  image,  calculated  from  Figs.  33  and 
34.  Figure  36  shows  the  average  vertical  profile  obtained  from  the  normalized  image,  showing 
again  a  ratio  of  about  0.6  between  the  FLIP  signal  below  and  above  the  centerline  in  the  FLIP 
image.  As  was  the  case  for  Run  14,  there  is  a  large  increase  in  signal  when  going  from  the  low- 
signal  region  to  the  high-signal  region.  From  Fig.  36,  this  ratio  is  calculated  to  be  22  (±2), 
consistent  with  the  value  found  for  Run  14. 

The  same  procedure  outlined  for  Run  14  was  followed  to  calculate  the  actual  position  of  the 
laser  sheet.  That  is,  a  set  of  points  was  selected  in  the  normalized  image  which  were  assumed  to 
define  the  intersection  of  the  laser  sheet  with  the  shock.  These  points  were  provided  as  input  to  a 
least-squares  fitting  program  to  optimize  the  laser  sheet  position.  Also  included  in  this  fit  was  the 
camera  pitch  angle,  5,  (Fig.  9),  because  it  appeared  that  the  camera  got  bumped  slightly  from  its 
original  alignment  position  (as  witnessed  by  a  vertical  shift  of  the«ymmetry  line  of  the  FLIP  image 
relative  to  the  pre-run  calibration  image  of  the  test  article). 

Fkguie  37  shows  two  sets  of  hand-picked  points  which  were  used  as  input  to  the  least-squares 
fitting  program.  Also  shown  are  the  optimized  intersections  of  the  laser  sheet  with  the  shock  for 
each  set  of  points.  Figure  38  shows  the  result  for  the  calculated  laser  sheet  positions  for  both  sets 
of  points.  Clearly,  the  calculated  positions  of  the  laser  sheet  depend  on  how  the  points  representing 
the  intersection  between  the  laser  sheet  and  the  shock  are  selected.  That  is,  how  is  the  position  of 
the  edges  defined  in  the  vertical  profile  of  Fig.  36?  However,  it  again  appears  that  the  laser  sheet 
moved  significantly  from  the  original  pre-run  position.  Moreover,  the  calculated  motion  of  the 
laser  sheet  for  Run  IS  is  almost  identical  to  the  calculated  motion  for  Shot  14. 

Figure  39  shows  the  result  of  remapping  the  normalized  PUF  image  to  a  median  plane 
through  the  test  article  (Fig.  32).  The  top  (a)  and  bonom  portions  (b)  ate  shown  separately.  Finally, 
Fig.  40  shows  an  enlargement  (calculated  separately)  of  the  detail  region  marked  in  Fig.  39b.  This  is 
the  region  that  contains  the  remapped  diagonal  “streak”  which  is  clearly  seen  in  the  left  bottom  comer 
of  the  raw  PLIF  image  from  Fig.  33  (a  less  prominent  streak  at  a  different  angle  is  also  visible  in  the 
top  left  region  of  the  image).  This  streak  is  likely  the  result  of  some  flow  disturbance  (perhaps  as  a 
result  of  a  temporary  d^osit  of  particulate  matter  on  the  nose  tip;  deposits  have  been  observed  after  a 
run).  As  a  result  of  the  remapping  transformation,  the  angle  of  this  streak  is  aligned  much  more  closely 
with  a  streamline  than  might  be  expected  from  the  roughly  4S-deg  angle  in  the  raw  FLIF  Image. 
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9.0  DISCUSSION  AND  SUMMARY 

The  number  of  PLIF  images  obtained  during  FY94  runs  of  the  Impulse  Facility  was  smaller 
than  anticipated.  One  consequence  of  this  was  that  no  images  were  obtained  for  identical 
operating  conditions  of  the  facility.  Nevertheless,  the  four  images  which  were  obtained  constitute 
a  sufficient  database  to  allow  a  detailed  comparison  with  CFD  predictions.  So  far,  such  a  detailed 
comparison  has  only  been  performed  for  the  PLIF  image  from  Run  16.  Although  incomplete  (for 
example,  the  effect  of  laser  attenuation  was  not  included  in  the  analysis  which  was  presented  in 
AIAA  paper  94-2621),  it  has  pointed  out  some  significant  discrepancies  between  the  measured 
and  calculated  data.  Quite  possibly,  these  discrepancies  point  at  inadequacies  of  the  CFD  model. 
However,  it  would  be  desirable  to  perform  further  experimental  testing  to  eliminate  the  possibility 
that  systematic  errors  are  introduced  in  the  measurement  process  and  the  PLIF  data  analysis. 
Specific  experiments  had,  in  fact,  been  planned,  but  had  to  be  forfeited  when  the  number  of 
Impulse  Facility  runs  was  curtailed  from  the  original  plan. 

In  addition  to  significantly  expanding  the  database  for  the  Impulse  Facility,  the  present  effort 
has  resulted  in  a  clear  understanding  of  the  process  which  is  required  to  acquire  quantitative  PLIF 
images  in  a  complex  facility.  This  includes  the  use  and  analysis  of  calibration  images  that  are 
required  for  proper  scaling  of  the  images.  For  example,  perfect  agreement  was  found  between 
measured  locations  of  the  bow  shock  determined  independently  by  PLIF  and  by  schlieien  imaging 
for  Run  16. 

No  such  agreement  could  be  established  for  Runs  14  and  15,  presumably  as  a  result  of 
motion  of  the  laser  sheet  during  the  Impulse  Facility  rans.  Clearly,  such  motion  is  undesirable 
and  should  be  avoided  in  future  tests  by  firmly  securing  the  laser  sheet  optics  (as  was  done 
following  Run  15).  However,  with  the  aid  of  the  PLIF  analysis  techniques  developed  during  this 
effort,  it  was  possible  to  reconstruct  the  likely  laser  sheet  positions  based  on  the  measured 
schlieren  shock  positions. 

Still  it  would  be  of  great  interest  to  repeat  the  setup  used  for  Runs  14  and  15  (laser  sheet 
passing  in  front  of  the  test  article),  because  the  behavior  of  the  PLIF  signal  is  dramatically  different 
from  that  in  Runs  16  and  19  (with  the  laser  sheet  along  a  meridian  plane  of  the  test  article).  By 
securing  the  laser  sheet  optics  securely  there  should  be  no  need  for  a  posttest  optimization  of  the 
laser  sheet  position.  An  analogous  geometry  of  interest  would  be  one  in  which  the  laser  sheet  is 
positioned  nearly  perpendicular  to  the  test  article,  so  that  the  intersection  of  the  laser  sheet  with  the 
shock  would  appear  as  a  crescent  arc. 

Having  successfully  developed  a  procedure  for  acquiring,  calibrating,  and  processing  of  PLIF 
images  in  the  Impulse  Facility,  the  emphasis  of  future  work  should  be  directed  in  particular  to 
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resolving  the  discrepancies  between  the  PLIF  images  and  the  computational  results.  The 
recommended  approach  for  this  is  to  select  a  particular  operating  condition  for  the  Impulse  Facility, 
and  acquire  a  set  of  PLIF  images  for  different  laser  sheet  geometries  and  different  excitation  wave¬ 
lengths  (perhaps  simultaneously  using  a  dual  PLIF  setup).  Along  with  this  experimental  work  there 
should  be  a  computational  effort  to  resolve  such  questions  as  to  why  the  predicted  nitric  oxide 
density  near  the  model  surface  might  be  too  low  or  why  the  predicted  temperature  in  this  region 
might  be  too  high.  The  combined  effort  should  lead  to  a  stronger  confldence  in  the  accuracy  of  the 
computational  results  and  might  lead  to  significant  improvements  in  modeling  the  complex  kinetics 
of  high  enthalpy  flows. 
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Dimensions  in  mm 

Figure  2.  Schematic  of  test  article. 


Side  View 


Top  View 


Side  View 


Top  View 


b.  Runs  16  and  19 

Figure  3.  Laser  sheet  geometries  relative  to  test  article. 
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Figure  5.  Result  of  digitizing  and  fitting  of  schlieren  shock  position  for 
Run  16. 


X,  mm 

Figure  6.  Fitted  schlieren  shock  positons  for  Runs  14, 15,  and  16. 
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Figure  7.  Image  of  calibration  target  used  for  Run  15. 
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Figure  8.  Image  of  calibration  target  used  fur  Run  16. 
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Figure  10.  Definition  of  camera  orientation  relative  to  test  article 
coordinate  system. 
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Figure  11.  Comparison  of  selected  and  fitted  marks  for  calibration  image 
from  Fig.  7. 
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Figure  12.  Overlay  of  fitted  test  article  features  on  calibration  image  from  Fig.  8. 
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Figure  13.  Examples  of  projections  of  remapped  image  pixel  onto 
original  image  grid. 
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Figure  14.  Flat-field  calibration  image  obtained  with  continuous  light  source 
and  long  gate. 
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Figure  15.  Typical  flat-field  calibration  image  obtained  with  a  short  laser  pulse 
and  a  short  gate. 
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Figure  17.  Run  16  PLIF  image  overlaid  with  projection  of  model  and  shock. 
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Divergence-Corrected  Horizontal  Coordinate 
Figure  18,  Average  horizontal  profile  of  free-stream  region  in  Fig.  16. 
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Figure  19.  Remapped,  energy-normalized  PLIF  image  for  Run  16. 
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Average  radial  PLIF  profiles  through  the  center  of  Fig.  19. 
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Figure  21.  Raw  PLIF  image  for  Run  19. 


Figure  22.  Remapped,  energy-normalized  PLIF  image  for  Run  19. 
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Figure  23.  Temperature  dependence  of  the  PLIF  signal  for  the  two  transitions 
used  in  this  work. 
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Figure  25.  Raw  PLIF  image  for  Run  14. 
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Horizontal  Pixel  Coordinate 

Figure  26.  Average  horizontal  profile  of  the  raw  image  from  Fig.  25. 
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Figure  29.  Expected  projections  of  the  intersection  of  the  laser  sheet  with 

the  shock,  as  seen  by  the  camera  for  Run  14,  for  various  distances 
between  laser  sheet  and  edge  of  the  model. 


Figure  30.  Overlay  of  selected  and  fitted  intersection  points  between 
laser  sheet  and  shock  on  energy-normalized  PLIF  image 
from  Fig.  27. 
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Schlieren  Shock 


Figure  31.  Fitter  laser  sheet  position  for  Run  14  shown  in  test 
article  coordinates. 
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Figure  32.  Remapped,  energy-normalized  PLIF  images  for  Run  14. 
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Figure  33.  Raw  PLIF  image  for  Run  15. 
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Figure  34.  Average  horizontal  profile  of  the  raw  image  from  Fig.  33. 


49 


AEDC-TR-94-14 


Figure  35.  Energy-normalized  PLIF  image  for  Run  15. 
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Figure  36.  Average  vertical  profile  of  the  PLIF  signal  in  Fig.  35. 
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Figure  37.  Overlay  of  selected  and  fitted  intersection  points  between  laser  sheet  and 
shock  on  energy-normalized  PLIF  image  from  Fig.  35. 
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Figure  38.  Fitted  laser  sheet  positions  for  Run  15  shown  in  test 
article  coordinates. 
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Figure  39.  Remapped,  energy-normalized  PLIF  images  for  Run  14. 
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Run  15:  Bottom  detail 
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Figure  40.  Detail  region  of  Fig.  39b. 
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Table  1:  Fitted  Coefficients  for  Schlieren  Shock  Positions 


*For  jr,  y  in  Eq.(2)  in  mm 
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NOMENCLATURE 

a,  b,  c  Fit  coefficients  for  schlieren  shock  position  from  Eq.  (2) 

Op  bj  Coefficients  giving  laser  sheet  position 

CCD  Charge-Coupled  Device  (camera) 

CFD  Computational  Fluid  Dynamics 

CFl  Computational  Row  Imaging 

d  Radial  distance  between  test  article  surface  and  laser  sheet 

D  Distance  between  camera  and  camera  reference  point  from  Fig.  10 

^x>  Camera  unit  vectors  defined  in  Eq.  (4) 

Enclosed  fraction  of  original  pixel  in  remapping  fiom  Eq.  (7) 

/.  j  Image  pixel  indices  in  remapping  process  expressed  by  Eq.  (7) 

ICCD  Intensified  Charge-Coupled  Device  (camera) 

^oHgineU  Original  image  pixel  intensity  from  Eq.  (7) 

^remapped  Remapped  image  pixel  intensity  from  Eq.  (7) 
p  Zoom  factor  for  camera  from  Eq.  (3) 

PLIF  Planar  Laser-Induced  Fluorescence 

R  Nose  radius  of  spherically  blunt  test  article 

T'^  Transformation  from  image  coordinates  to  test  article  coordinates 

w  Vector  defined  in  Eq.  (5) 

X,  y,  z  Coordinates  associated  with  test  article  reference  frame  (Fig.  10) 
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Xp,  yp  Image  pixel  coordinates 

Xo>  ych  ^  Coordinates  of  point  at  which  camera  is  pointed 

a  Half-angle  of  conical  surface  of  test  article 

Y  Off-peipendicular  camera  viewing  angle  from  Fig.  10 

8  Camera  pitch  angle  from  Fig.  10 

^  Fit  coefficient  for  schlieren  shock  position  from  Eq.(2) 

(p  Camera  roll  angle 
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